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Recently we have shown that the centrality dependence of the multiplicity of different hadron 
species observed in RHIC and SPS experiments can be well understood in a simple model, dubbed 
core-corona model. There it is assumed that those incoming nucleons which scatter only once 
produce hadrons as in pp collisions whereas those which scatter more often form an equilibrated 
source which decays according to phase space. In this article we show that also kinematical variables 
like V2 / epart{Npart) as well as v\/ epart[Npart) and < p^j,[Npart) > of identified particles are well 
described in this model. The correlation of < > between peripheral heavy ion collisions and 
pp collisions for different hadrons, reproduced in this model, questions whether hydrodynamical 
calculations are the proper tool to describe non-central heavy ion collision. The model explains as 
well the centrality dependence of V2 / epart of charged particles, considered up to now as an observable 
which allows to determine the viscosity of the quark gluon plasma. The observed dependence of 
V2/ (-part{Npart) On the particle species is a simple consequence of the different ratios of core to 
corona particles. 
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Simulations of heavy ion collisions with advanced event 
generators like EPOS [? ], which reproduce a multitude 
of experimental observables, have revealed that nucleons 
at the surface of the reaction zone (called corona parti- 
cles) have only few collisions and do not come to statis- 
tical equilibrium with the more central (core) particles 
which form an equilibrated system. To study the conse- 
quences of this observation we have developed a simple 
model by defining corona particles as those nucleons 
which have only one initial collision whereas the others 
are considered as core particles, fcore is the fraction of 
core nucleons which depends on the centrality, the sys- 
tem size and (weakly) on the beam energy. In this sim- 
ple model we could show that, independent of the system 
size, the centrality dependence of the multiplicity of all 
hadrons from SPS to RHIC energies can quantitatively 
be described by: 
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where fcore, shown in Fig. [TJ has been calculated in a 
Glauber model and i refers to the hadron species. For our 
calculation we fix M^^^^ by applying eq. [1] to the most 
central AuAu or PbPb data point. Ml^^^^^ is given as 
half of the multiplicity measured in pp collisions. Once 
these parameters are fixed, the centrality dependence of 
is determined by eq. [TJ Especially the centrality de- 
pendence of the lighter CuCu system follows then with- 
out any further input. Certainly this is a very simple 
model with no free parameter (besides the Glauber cal- 
culation of fcore) but the present experimental error bars 
of the quantities which we analyze give not sufficient in- 
formation to refine the model. 



AuAu/CuCu 200 AGeV 



1.0 
0.9 
0.8 
0.7 

§ 0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



n O 

^00 



O 



o 



o 



AuAu 
CuCu 



50 100 150 200 250 300 350 400 

Npart 



FIG. 1: fcore the fraction of core nucleons as a function of 
the participant number Npart for AuAu and CuCu collisions 
at Vs = 2mGeV. 



One may ask whether the core - corona model can also 
describe the centrality dependence of other observables, 
like < pt > (Npart) or V2/ tpart{Npart) of charged parti- 
cles or identified hadrons. Especially the centrality de- 
pendence of V2 / epart{Npart) has recently created a lot of 
theoretical activities. Initially the azimuthal distribution 
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da/dcj) oc (1 + 2vi ■ cos((/)) + 2v2 ■ cos(20)) of light quarks 
and gluons is isotropic and the anisotropy develops dur- 
ing the expansion as an image of the initial eccentricity 
£part in coordinate space. This initial eccentricity has 
been calculated in a Glauber model or in a Color Glass 
Condensate approach. Taking this eccentricity as input, 
calculations based on ideal hydrodynamics reproduce the 
V2 Impart observed in central AuAu collisions at RHIC. 
They predict, however, a centrality independent V2/epart 
value and fail therefore to describe the measured central- 
ity dependence of Vil^part which decreases by a factor of 
two from central to peripheral reactions. 

Drescher et al. Q argue that this centrality depen- 
dence is a sign that the parton cross section in the quark 
gluon plasma, created presumably at RHIC energies, is 
too small in order to justify the application of ideal hy- 
drodynamics. It has to be replaced by viscous hydro- 
dynamics. They could show that the impact parameter 
dependence of vijf-part can be related to the Knudsen 
number and therefore to a finite viscosity. Later Luzum 
et al. [3] showed in viscous hydrodynamical calculations 
that different viscosities yield a different centrality de- 
pendence of V2 Impart ^ud of < >■ Performing calcu- 
lations with different values for the viscosity they could 
identify the viscosity which fits best the 200 AGeV AuAu 
data and propose that this value is characteristic for the 
quark gluon plasma. More recently. Song et al. 5] have 
extended the calculation to the CuCu system. It is there- 
fore challenging to explore the predictions of the core - 
corona model for dynamical observables and to compare 
the results to data. 

The centrality dependence of < > at midrapidity 
predicted by the core - corona model is given by 

<Pt> (^part) (2) 
~ fcore' ^ Pt >corc +(1 ^ fcore)' < Pt >corona ■ 

where i characterizes the particle species i. < pip >corona 
is the value of < > measured in pp collisions at midra- 
pidity whereas < pip >coro is the value of < > mea- 
sured in central heavy ion reactions at midrapidity. /^o^e 
is the centrality dependent fraction of core particles for 
the particles species i, given by 

n _ J core -'"core ('o^ 

~ f . A,fi _L n _ f 1 . ■ ' 

J core ^"corc ^ \^ J core J -'"corona 

If the system consists of particles of one species, final 
state interactions change neither < px > nor V2/epart 
substantially. The former is due to (transverse) en- 
ergy conservation, the latter due to the small momen- 
tum transfer in final state interactions as compared to 
the momentum of the particles. It is also not important 
whether some corona particles are absorbed by the core. 
This increases the core fraction but simultaneously low- 
ers < pt >core- Both cffccts Compensate as long as the 
total transverse energy is conserved. 

This observation is not anymore true if the systems 
contains different particles species i. As we will see, the 



dependence of < > on Npart is quite different for dif- 
ferent hadrons. Interactions transfer momentum between 
the different particles species and hence change < pip > 
in the direction toward a thermal equilibrium. Therefore 
we do not expect that < pip > follows eql2] if final state 
interactions between the hadrons become important. 

Consequently, if eq. [5] is in agreement with data, we 
can draw conclusions on the reaction scenario. 

The NA49 collaboration Q has recently analyzed the 
transverse mass spectrum of different hadrons in the 
framework of our model and has found good agreement 
with data. Therefore it is meaningful to continue here 
this analysis towards higher energies. 
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FIG. 2: < > as a function of the centrality for different 
particle species i. The STAR data fMH^ are presented by 
symbols, the calculations, eql2] by lines. Full (open) symbols 
refer to CuCu (AuAu). The < pr > values of protons (cj)) are 
multiplied by a factor of 1.75(1.4). 

In fig. [5] we compare our predictions with published 
0-il and preliminary [13, [U STAR data, in fig. [3] with 
PHENIX data In fig. [2] the theoretical and experi- 
mental values of protons {(j)) are multiplied by a factor of 
1.75(1.4) to separate the curves. The STAR proton data 
show (in contradistinction to the proton data of PHENIX 
and to those of all other hadrons) a strong increase of 
< Pt > for Npart > 100. We have therefore multiplied 
the experimental < pr >core by 0.9. 

We see that < plp{Npart) >, calculated with help of eq. 
m agrees quantitatively with the data for the different 
hadrons. Also the experimentally observed different de- 
pendencies of < plp{Npart) > for different particle species 
i is well reproduced. We note in passing that replacing 
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TABLE I: A summary of mid-rapidity (pt) in GeV/c for 
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FIG. 3: < > as a function of the centrality for different 
particle species i. The PHENIX data [l^ are presented by 
symbols, the calculations by lines. 



fcore by fcore in eq. [2] changes the curve only slightly. Its 
form is dominated by the difference between < pT >core 
and < Pt >corona- TabU summarizes the experimental 

< Pt > corona and < p^rp >corc which enter our calcula- 
tions and which are taken from pp and central heavy ion 
collisions. PHENIX has measured very peripheral heavy 
ion collisions. The < pT > value there coincide within 
the error bars with the STAR pp data. We note that 

< Pt > corona / < Pt >corc varies from .81 for pions to 
.62 for protons. 

By construction the core - corona model correlates the 

< > values of peripheral heavy ion reactions with 
those observed in pp. This correlation is shown in figlH 
It displays < p^j, >peri / < Pt >core as a function of 

< Pt > corona / < Pt >core, usiug cqlSl for different 
hadron species. The experiential data show this correla- 
tion as well, it is even stronger there. Thus also in exper- 
iment the value of < pip > in peripheral reactions is small 
as compared to that in central collisions for those hadrons 
for which < pip > measured in pp is small as compared 
to < > measured in central heavy ion collisions. This 
correlation is far from being trivial. If, as often assumed, 
the heavy ion collision can be described by hydrodynam- 
ical calculations we do not expect any correlation be- 
tween peripheral heavy ion data and pp data. Therefore 
the correlations observed in fig. |4] support the validity 
of the core - corona approach and questions whether hy- 
drodynamical calculations are the proper tool to describe 
non-central heavy ion reactions. 
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FIG. 4: The ratio of < pT > in peripheral to central reactions 
as compared to that of pp to central reactions. The data are 
from [7 h1i1 . [i3 |. Because the calculation uses the experimental 
data, the error bars of model predictions are of about the same 
size as that of the experimental points. 



Corona particles do not form a thermal equilibrium 
with the core particles but decay like pp collisions and 
hence isotropically. Only the core particles feel the eccen- 
tricity of the overlap region. Because M^j^^^^'^'^ particles 
M^orolf V2/epart is cxpected to be oc Ure and 

therefore the expected centrality dependence of vije-part 
is that of fcore- This allows for some immediate pre- 
dictions. Because fcore is similar in central CuCu and 
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AuAu collisions (Fig. [Ij we expect that central collisions 
of AuAu and CuCu show a similar v^jf-yart despite of 
the large difference of Npart- In peripheral CuCu and 
AuAu collision V2/epart should be very similar for the 
same Npart- In order to compare our model with exper- 
iment we have to determine {v2 / tpart)^^'''^" ^ i-e. v^jtpart 
for /core = 1. We take ivij tT>artY^'^''° as 24%, 19.7% 
and 12.6% for reactions at [fs of 200 AGeV, 62 AGeV 
and 17.2 AGeV, respectively. In figl5]we display the ex- 
perimental results from PHOBOS (full symbols), from 
STAR (open symbols) and from NA49 (open squares) 
in comparison with our predictions, presented as lines. 
We see a quite nice agreement. As in our model, for a 
given number of Npart the experimental values of vi / tpart 
for CuCu are higher than those for AuAu. Hydrody- 
namical calculations predict the opposite. This allows 
to conclude on the underlying mechanism once the error 
bars of the present preliminary data are further reduced. 
The situation becomes much more clear when we divide 
y2/tpart by /core- Our modcl predicts a constant. Fig. [S] 
presents the PHOBOS data for the 4 different systems, 
AuAu and CuCu at 200 and 64 AGeV . The experimen- 
tal V2 / impart f core) IS iudccd independent of Npart which 
means that the centrality dependence of V2/{epart) agrees 
with that of fcore- This does not exclude, of course, 
that in hydrodynamical calculations a viscosity can be 
obtained by fitting the centrality dependence of the ex- 
perimental data. In our approach we do, however, not 
need such a fit parameter. It is astonishing that even for 
the most peripheral reactions, where the number of core 
particles is of the order of 10 and therefore the assump- 
tion of a hydrodynamical elliptical flow is questionable, 
the predictions of core - corona model still hold. 

In our approach V2ltpart is a property of the plasma 
particles and therefore one may expect the same central- 
ity dependence of V2/€part for all hadrons (although we 
cannot exclude that the hadronization itself creates dif- 
ferences for different types of hadrons). This is, how- 
ever, not true. The fraction of particles which come 
from the core depends, according to eq. [1] on the par- 
ticle species. For all charged particles and tt, for which 
Mcore ~ Afcorona, this dependence is weak jl^ (and there- 
fore we have divided V2 by fcore) but for multi-strange 
hadrons, for which Mcoxc » A/corona, this is an impor- 
tant effect. We expect 
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nary pt integrated Kp, A, S data as measured by the 
STAR collaboration [13]- These data are compared 
with the model, eq. 2) Model and experiment agree 
within error bars and the trend of the data is well 
reproduced. We note that the pt integrated value of V2 
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FIG. 5: Top: Npart dependence of 112 /tpart for charged 
hadrons in AuAu and CuCu at 200 AGeV and 62 AGeV as 
measured by PHOBOS [3 and STAR [H, [iTl] as well as at 
17.2 AGeV [H] as measured by NA49 in comparison with 
theory. Bottom: Zoom of the region of small Npart- 



does not show a scaling with the number of entrained 
quarks Uq which has been observed in the differential 
spectra for (rriT — mo)/nq < IGeV , with ttit (mo) being 

the transverse (rest) mass. The relative ratio of — ^ 
between baryons and mesons depend on centrality and 
on the particle species and is not 3/2 corresponding 
to the ratio of entrained constituent quarks. However, 
before firm conclusions can be drawn more data have to 
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FIG. 8: EPOS calculation compared with data. The top panel 
shows 112 as a function of pr for Tr**" measured in semi-central 
reactions. The bottom panel shows the pt integrated value 
of V2 as a function of the number of participants. Data are 
from ref. p^j . 
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Pt integrated Ko,A,S. data as measured by the STAR 
collaboration [l7| in comparison with the prediction of the 
core-corona model, eq|4l presented as lines. 



be available. Especially data for protons for which /^ore 
and f corona differ much less than for the multi-strange 
hadrons would be of help. Such an analysis will not be 
easy as can be seen from fig. [8] which shows V2 of tt's 
as a function of pr for semi-central reactions (20%-60% 
centrality) (top panel) and that of charged particles as 
a function of Npart (bottom panel) . We compare results 
of EPOS calculations with the experimental data. The 
top panel shows that the V2{pt) spectra is qualitatively 
described up to pr ~ IGeV/c, the range where almost 
all particles are located. Nevertheless, after integration 
over Pt, bottom panel, we see differences of V2 of the 
order of 10% between theory and experiment, although 
V2 of tt's and of charged particles should be very similar. 
This is due to the extrapolation of the experimentally 
measured curve toward px — >■ 0, necessary to determine 
a Pt integrated V2 value. 

Summarizing, the centrality dependence of all observ- 
ables for identified particles, their multiplicity, their av- 
erage transverse momentum and their elliptical flow, as 
well as the elliptical flow of charged hadrons can be quan- 
titatively described from SPS to RHIC energies in the 
very simple core - corona model. It would be interesting 
to compare these results in detail with the predictions 
of (viscous) hydrodynamical calculations, the alternative 
approach to describe non central heavy ion collisions. 
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The model predicts even the experimentally observed 
correlations between peripheral heavy ion and pp colli- 
sions. Such a correlation is alien to any hydrodynamical 
approach. Also the larger W2/epart value for CuCu as 
compared to AuAu for the same number of Npart ques- 
tions whether the centrality dependence of V2/€part can 
be related to the viscosity of the plasma. 

The success of the core - corona model suggests the 
following reaction scenario: The experimental data are 
a superposition of two contributions: a corona contribu- 
tion which is given by pp physics and a core contribution 
which agrees with the assumption that the core forms an 
equilibrated source. The data are compatible with the as- 
sumption that the properties of this source are indepen- 



dent of centrality. A centrality dependence of observables 
appears because the relative fraction of core and corona 
contribution varies with centrality. Core particles are 
in statistical equilibrium. Their interaction with corona 
particles as well as the interaction among corona particles 
has to be unimportant otherwise momentum would be 
transferred between the different species and eq. [2] would 
not be valid anymore. This suggest also a very limited 
number of collisional interactions between hadrons from 
core and corona. 
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